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Genetic susceptibility to carbamazepine-induced cutaneous 
adverse drug reactions 

Shuen-lu Hung a \ Wen-Hung Chung b,c ' 1 , Shiou-Hwa Jee d , Wen-Chieh Chen e , 
Yun-Ting Chang f , Woan-Ruoh Lee 9 , Shu-Ling Hu h , Meng-Tse Wu', 
Gwo-Shing Chen j , Tak-Wah Wong k , Pa-Fan Hsiao', Wei-Hsuan Chen 3 , 
Han-Yu Shih a , Wu-Hsiang Fang 3 , Chun-Yu Wei 3 , Yi-Hui Lou 3 , Yau-Li Huang b , 
Juei-Jueng Lin m and Yuan-Tsong Chen a,n 



The anticonvulsant carbamazepine (CBZ) frequently 
causes cutaneous adverse drug reactions (cADRs), 
including maculopapular eruption (MPE), hypersensitivity 
syndrome (HSS), Stevens-Johnson syndrome (SJS) and 
toxic epidermal necrolysis (TEN). We reported that 
SJS/TEN caused by CBZ is strongly associated with the 
HLA-B*1502 gene in Han Chinese. Here, we extended our 
genetic study to different types of CBZ-cADRs (91 patients, 
including 60 patients with SJS/TEN, 13 patients with 
hypersensitivity syndrome and 18 with maculopapular 
exanthema versus 1 44 tolerant controls). We used 
MALDI-TOF mass spectrometry to screen the genetic 
association of 278 single nucleotide polymorphisms 
(SNPs), which cover the major histocompatibility complex 
(MHC) region, tumor necrosis factor-alpha, heat shock 
protein and CBZ-metabolic enzymes, including CYP3A4, 
2B6, 2C8, 2C9, 1A2 and epoxide hydrolase 1. In addition, 
we genotyped 20 microsatellites in the MHC region and 
performed HLA-typing to construct the recombinant map. 
We narrowed the susceptibility locus for CBZ-SJS/TEN to 
within 86 kb flanking the HLA-B gene on the extended 
B*1502 haplotype, and confirmed the association of 
B*1502 with SJS/TEN [P c =1.6 x 10" 41 , odds ratio 
(OR) = 1357; 95% confidence interval (CI) = 193.4-8838.3], 
By contrast to CBZ-SJS/TEN, HLA-B* 1502 association was 
not observed in the MPE or HSS groups: MPE was 
associated with SNPs in the HLA-E region and a nearby 
allele, HLA-A*3101 (P c =2.2 x 10" 3 , OR = 17.5; 95% 
CI =4.6-66.5), and HSS with SNPs in the motiiin gene 
(P c =0.0064, OR = 7.11; 95% Cl = 3.1-16.5) located terminal 
to the MHC class II genes. No SNPs in genes involved in 
CBZ metabolism were associated with CBZ-induced 
cADRs. Our data suggest that HLA-B*1502 could 
contribute to the pathogenesis of CBZ-SJS/TEN, and that 
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Introduction 

Although carbamazepine (CBZ) is a commonly prescribed 
first-line anticonvulsant, it is also a relatively common 
cause of cutaneous adverse drug reactions (cADRs) [1]. 
The frequency of CBZ-induced cADRs is between 1 in 

'Both authors contributed equally to this work. 
1744-6872 © 2006 Lippincott Williams & Wilkins 



1000 and 1 in 10000 new exposures to the drug in 
Caucasians [2,3], The cADRs range from mild maculo- 
papular eruption (MPE), with increasing severity, to 
hypersensitivity syndrome (HSS), Stevens-Johnson syn- 
drome (SJS) and toxic epidermal necrolysis (TEN). The 
mortality rate of the latter can be as high as 40% [4]. 
MPE is characterized by cutaneous fine pink macules and 
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papules, lesions which usually fade within 1-2 weeks 
following cessation of drug treatment. HSS is character- 
ized by multi-organ involvement (e.g. hepatitis and 
nephritis) accompanied by systemic manifestations (e.g. 
fever, arthragias, eosinophilia and lymphadenopathy) in 
addition to skin rashes [4). Skin manifestations of HSS 
may vary from MPE to exfoliative dermatitis [5]. HSS is 
also called drug rash with eosinophilia and systemic 
symptoms. SJS and TEN are characterized by a rapidly 
developing blistering exanthema of purpuric macules and 
target-like lesions accompanied by mucosal involvement 
and skin detachment [6]. According to the clinical 
morphology, SJS/TEN belong to the group of bullous 
cADRs, whereas MPE and HSS are non-bullous reactions. 

The mechanism by which CBZ causes cADRs is not well 
understood. Potential defects in the enzymes responsible 
for bioactivation and detoxification of CBZ have been 
proposed [1]. Carbamazepine is bioactivated by hepatic 
cytochrome P450 enzymes, mainly CYP3A4, CYP2B6 and 
CYP2C8, which generate various potentially reactive 
metabolites, such as CBZ 10,1 1-epoxide, 3-hydroxy 
CBZ, 2-hydroxy CBZ, and CBZ 2,3-epoxide [7,8). Most 
of the reactive epoxides are detoxified to non-toxic 
dihydrodiols by liver microsomal epoxide hydrolase 1 
(EPHX1) or to glutathione conjugates by glutathione 
transferase [9,10]. Previous attempts to identify a genetic 
defect altering the structure or function of epoxide 
hydrolase 1 in individuals susceptible to CBZ-induced 
cADRs, however, were not successful [11,12]. 

It has also been proposed that immune reactions are 
involved in the cADRs triggered by CBZ because 
infiltrating inflammatory cells can be detected in the 
skin lesions, and rechallenging with the same drug 
typically shortens the incubation period, resulting in 
more severe manifestations [13,14]. CD4 + T cells are 
the major cell type found in the skin lesions of MPE 
and HSS [15], whereas CD8 + T-cell-mediated 
cytotoxic responses appear to be the major event in 
SJS/TEN [16], There is also evidence supporting the 
view that cADRs involve major histocompatibility com- 
plex (MHC)-dependent presentation of its metabolites 
forTcell activation [16,17]. Naisbitt et al. [15] reported 
that CBZ might be presented by MHC class II expressed 
on the surface of antigen-presenting cells to the T-cell 
receptor of CD4 + T cells of patients with CBZ-induced 
HSS [15]. 

We recently reported that, in Han Chinese, the HLA- 
B* 1502 gene is strongly associated with CBZ-induced SJS 
and TEN [18]. In the present study, we sought to extend 
our genetic association study to other types of cADRs 
induced by CBZ using a comprehensive genetic screen of 
markers in the MHC region, as well as genes encoding 
CBZ metabolic enzymes. In addition, we performed fine 



mapping to further investigate the susceptibility locus for 
CBZ-induced SJS/TEN. 

Materials and methods 
Study population 

From 1997 to 2004, 91 individuals who fulfilled diagnostic 
criteria for CBZ-induced MPE, HSS, SJS or TEN agreed 
to participate in the study. Eighty-eight of the 91 patients 
were recruited from Chang Gung Memorial Hospital 
Health System, Veterans General Hospital, National 
Taiwan University Hospital, Cathay General Hospital 
and Chushang Show-Chwan Hospital, Taiwan. We 
recruited patients primarily from inpatients who had 
more severe diseases. The remaining three patients were 
referred by physicians from Hong Kong (two patients) 
and the USA (one patient). The latter three all had CBZ- 
SJS. All patients (including the three from Hong Kong 
and the USA) and controls (see below) were Han 
Chinese or Chinese descendants. 

The 91 patients included 44 patients with CBZ-SJS/TEN 
who we previously reported [18], and 16 additional 
patients with SJS/TEN, as well as 13 patients with 
hypersensitivity syndrome and 18 with maculopapular 
exanthema. All patients were assessed by two dermatol- 
ogists who reviewed photographs, pathological slides and 
medical records. Diagnostic criteria for SJS/TEN were 
based on the clinical morphology defined by Roujeau [6]. 
We defined 'SJS' as a skin detachment of less than 10% of 
total body-surface area, 'overlap of SJS and TEN' as skin 
detachment of 10-30%, and 'TEN' as skin detachment 
greater than 30%. The criteria for HSS in this study were 
skin rash, plus two of the following symptoms: fever, 
lymphadenopathy and haematologic abnormalities (e.g. 
eosinophilia, atypical lymphocytosis) with involvement of 
at least one internal organ (e.g. hepatitis, pneumonitis, 
myocarditis, pericarditis, nephritis) [2,5,19], SJS and 
TEN are bullous cADRs and are considered to be variants 
of the same disease. Therefore, they were analysed 
together as a group. MPE and HSS are non-bullous and it 
is less clear whether they are from a single disease 
spectrum. They were therefore analysed as a group and 
also analysed separately as two entities. In all enrolled 
cases, CBZ was regarded as the offending drug if the 
onset of cADRs symptoms occurred within the first 2 
months of exposure and the symptoms resolved upon 
withdrawal of the drug. Patients with an absence of 
symptoms after re-exposure to CBZ were excluded. 

The control group was the 144 consecutive patients who 
received CBZ for at least 3 months without evidence of 
adverse reactions. These tolerant patients were recruited 
from the Neurology Clinic of the same regional hospitals 
where cADRs patients were recruited. In addition, 93 
healthy subjects were randomly selected from a biobank 
under a nationwide population study, in which 3312 Han 
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Chinese descendants were recruited based on the 
geographical distribution across Taiwan. There was no 
self-report of adverse events in any of these 93 subjects. 
The study was approved by the institutional review 
board, and informed consent was obtained from all of the 
participants. 

DNA isolation and genotyping 

Genomic DNA was isolated using the PUREGENE DNA 
purification system (Gentra Systems, Minneapolis, Min- 
nesota, USA). 

Short tandem repeat polymorphism (STRP) genotyping 

Twenty highly polymorphic microsatellite markers 
located in the MHG region were selected from the 
NCBI database (i.e. D6S258, D6S2972, D6S510, 
D6S265, D6S388, D6S2814, HLAC-CA, HLABC-CA, 
MIB, MICA, TNFd, BAT2-CA, D6S273, D6S1615, 
DQCAR, G51152, D6S2414, D6S1867, D6S1560 and 
D6S1583). The average heterozygosity of markers was 
0.702 with an estimated 230 kb of spacing. Primers were 
designed based on oligonucleotide sequences reported 
within the database. Polymerase chain reaction (PGR) for 
genotyping was performed in a 5-ul volume containing 
lOng of genomic DNA and 0.33 um of each primer by 
using GeneAmp 9700 thermocyclers (Applied Biosys- 
tems, Foster City, California, USA). Up to six products of 
appropriate size and carrying a fluorescent label were 
pooled for capillary gel electrophoresis. The size of 
polymorphic amplicons was determined by ABI 3730 
DNA sequencer (Applied Biosystems) using the LIZ500 
size standard as an internal size standard. Allele sizing 
was calculated using the GENMAPPER program version 
3.0 (Applied Biosystems). Allele calling and binning were 
performed using the SAS program (SAS Institute, Cary, 
North Carolina, USA). Three CEPH control individuals 
(1331-01, 1331-02, 1347-2) and H 2 0 were included in all 
genotyping experiments for quality control purposes. 

Single nucleotide polymorphisms (SNP) genotyping on 
MHC region and metabolic enzymes for CBZ 

A total of 379 SNPs were selected from the SNP database 
(dbSNP: build 123) of the National Center for Biotech- 
nology Information for genotyping [20]. These included 
220 SNPs from 4Mb of the MHC region on chromosome 
6p21.3 and 159 SNPs selected from genes encoding drug 
metabolizing enzymes. The 220 SNPs selected in the 
MHC region included 201 SNPs reported by Walsh et al. 
[21], and additionally included SNPs of tumor necrosis 

factor (TNF) -alpha and the heat shock proteins. We 
especially included rsl800629, TNF-308A allele, which is 
reported to be associated with CBZ-HSS [22], and 
rs2227956 of HspAlL (heat shock 70 kDa proteinl-like), 
which is reported to be associated with abacavir 
hypersensitivity [23]. The average space of selected 
SNPs in the MHC region was approximately 20 kb. 
Another 159 SNP set contained genes potentially 



involved in bioactivation or detoxification of CBZ 
metabolism, including CYP3A4, CYP2B6, CYP2C8, 
CYP2C9, CYP1A2 and EPHX1. 

SNP genotyping was performed using high-throughput 
MALDI-TOF mass spectrometry. Briefly, primers and 
probes were designed using SpectroDESIGNER software 
(Sequenom, San Diego, California, USA). Multiplex 
PCRs were performed, and unincorporated dNTPs were 
dephosphorylated using shrimp alkaline phosphatase 
(Hoffman-LaRoche, Basel, Switzerland) followed by 
primer extension. The purified primer extension reaction 
was spotted onto a 384-element silicon chip (Spectro- 
CHIP, Sequenom), analysed in the Bruker Biflex III 
MALDI-TOF SpectroREADER mass spectrometer (Se- 
quenom) and the resulting spectra processed with 
SpectroTYPER (Sequenom). 

Of the 379 SNPs selected from the database, 101 SNPs 
were excluded from further analysis due to the non- 
polymorphism, low successful rate (less than 0.7), and/or 
departure from Hardy-Weinberg equilibrium (P < 0.001) 
in our control group. The remaining 278 SNPs, including 
88 SNPs in metabolic pathway for CBZ and 190 SNPs in 
the MHC region, were used for the statistical analyses 
and the P-values of association were adjusted by using 
Bonferroni correction. 

HLA genotyping 

HLA alleles A, B, C and DRB1 were determined by 
sequence-specific oligonucleotide reverse lineblots (DY- 
NAL Biotech Ltd, Bromborough, UK) [24]. Potential 
ambiguities were resolved by sequencing-based typing 
[25]. The primers used for PCR amplification of DNA 
fragments of HLA-A, B and C genes were: (i) Binl-TA- 
M 13F, Bin 1 -CG-M 1 3F and Bin3-M 13R for HLA-B alleles; 
(ii) Ainl-A-M13F, Ainl-G-M13F, Ainl-T-M13Fand Ain3- 
62-M13R for HLA-A alleles; and (iii) 5Clnl-61 and 
3BCln3-12 for HLA-C alleles [26]. The DRB1 alleles were 
separated into eight groups: DR1, DR2, DR3/11/6, DR4, 
DR7, DR8/12, DR9, DR10 for sequencing-based typing 
as described [27]. Genotyping for MHC class I chain- 
related gene A (MICA) was carried out by DNA 
sequencing on its exons 2, 3, 4 and 5 [25]. The 
primers used for PCR amplification of the DNA fragment 
of the MICA gene were: forward primer: 5'- 
CGTTCTTGTCCCTTTGCCCGTGTGC-3' and re- 
verse primer: 5'-GATGCTGCCCCATTCCCTTCCCAA-3'. 
The sequence data were analysed by SeqScape v2 
(Applied Biosystems). 

Statistical analysis 

Comparisons of allele or genotype frequencies between 
groups were performed using Fisher's exact tests. All 
/^-values were two-tailed. /"<0.05 was considered to be 
statistically significant. An allelic association screen was 
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carried out by the Cochran-Armitage Trend test for each 
SNP or STRP marker [28]. To test the association of 
haplotype frequencies, the composite haplotype method 
and haplotype trend regression of Helix Tree software 
version 3.0.0 were used. The linkage disequilibrium 
between two loci was analysed by using Expectation/ 
Maximization (EM) method provided by Helix Tree 
software version 3.0.0 (Golden Helix Inc., Bozeman, 
Minnesota, USA). To achieve sufficient power to identify 
loci associated with different clinical manifestations, the 
corrected P (P c ) values were adjusted by using Bonferro- 
ni's correction for multiple comparisons (278 for SNP 
assays, 17 for HLA-A, 40 for HLA-B, 19 for HLA-C, 30 for 
HLA-DRB1, 25 for MICA, and 20 for HLABC-CA). Odds 
ratios (ORs) were calculated using Haldane's modifica- 
tion, which adds 0.5 to all cells to accommodate possible 
zero counts [29]. 

Results 

Characteristics of patients and controls 

Based on the clinicopathologic features of cADRs, our 
patients could be divided into bullous and non-bullous 
groups: (i) bullous: CBZ-SJS/TEN (» - 60; 54 with SJS, 
five with overlapping SJS/TEN, one with TEN) and (ii) 
non-bullous: CBZ-MPE/HSS (» = 31; 18 with MPE, 13 
with HSS). The non-bullous group was further separated 
into two subgroups (MPE and HSS) for data analysis. 
Clinical manifestations and demographic variables of the 
91 patients and 144 tolerant controls are summarized in 
Table 1. The onset of symptoms for all patients with 
cADRs occurred within the first 2 months of CBZ 
exposure. The mean duration of CBZ exposure was longer 
in HSS patients than in SJS/TEN. Three patients with 
SJS/TEN had a second attack within 2 days of re- 
exposure; one of them developed TEN and died during 
the second attack. All 60 patients with SJS/TEN had 
widespread purpuric rash with blisters and mucosal 
involvement. In addition to maculopapular exanthema, 



all 13 patients with HSS had symptoms of multi-organ or 
systemic involvement. 

One hundred and forty-four patients who had been on 
CBZ for at least 3 months (mean 75.03 months) with no 
self-report of adverse events were used as tolerant 
controls. Their demographic variables are shown in Table 
1. Note that the CBZ-tolerant controls received higher 
doses of CBZ (mean dosage 793.8 mg/day) and yet no 
adverse drug reactions were observed. 

Fine mapping of the major histocompatibility 
complex region for genetic susceptibility to 
CBZ-SJS/TEN 

To confirm our previous observation of the genetic 
association between HLA-B*1502 and CBZ-SJS/TEN, 
and to further define the susceptibility region, we 
genotyped 20 STRP markers and 190 SNPs in the 
MHC region. Nine STRP markers located between 
D6S2814 and D6S2414 (physical map 30.8-32.9 Mb, 
chromosome 6) showed significant association 
CP-values < 10~ 6 ; Fig. 1). In particular, HLABC-CA near 
HLA-B showed the strongest association (P = 3.4 x 10" 19 ; 
Fig. 1), followed by D6S1615, a marker closed to HLA- 
DRB1 locus. 

Thirty-seven SNPs showed /'-values of less than 0.001 
(Fig. 2a). Among them, seven SNPs located between 
HLA-DRA and HLA-C (physical position 31.3-32.5 Mb, 
chromosome 6) showed strong associations (P< 10" l0 ) 
(Fig. 2a). The three most important SNPs were 
rs3130690, rs2848716 and rs750332, located near HLA- 
B, MICA and BAT3 genes (Table 2). The rs3130690, an 
intergenic SNP with 36 kb telomeric to the HLA-B locus, 
demonstrated the strongest association with CBZ-in- 
duced SJS/TEN (P c = 1.29 x 10" w ) (Table 2). The T 
allele of the rs3130690 SNP was present in 95% (57/60) 
of CBZ-SJS/TEN patients, but only in 6.9% (10/144) of 



Table 1 Demographic variables, dosage and duration of carbamazepine (CBZ) exposure, and clinical characteristics in CBZ-induced 
cutaneous adverse reactions and CBZ-tolerant controls 





Bullous cADRs 


Non-bullo 


us cADRs (n=31) 


Tolerant controls (n=144) 


SJS/TEN (n = 60) 


HSS (n = 13) 


MPE (n = 18) 


Sex, n (%) 
Male 










33 (55) 


9 (69.2) 


8 (44.4) 


75 (52.1) 




27 (45) 


4 (30.8) 


10 (55.6) 


69 (47.9) 


Age (years), mean (range) 


43.4 (5-80) 


51.5 (6-83) 


45.9 (7-84) 


35.7 (5-79) 


CBZ exposure, mean (range) 








793.8 (100-1500) 


Dosage (mg/day) 


332.6 (1 00-600) 


420 (300-600) 


423.1 (200-800) 


Duration 


15.1 (2-49) days 


32.7 (15-54) days 


22.4 (7-55) days 


75.03 (3-287.5) months 


Cutaneous features, n (%) 










Blister or epidermal detachment 




O(O) 


0 (0) 


Not observed 


Mucosal erosions 


60 (1 00) 


3(23) 




Not observed 


General and laboratory, n (%) 










High fever (>38.5°C) 


41 (68.3) 


12(92.3) 


11 (61.1) 


Not observed 


Eosinophil count > 1 000/mm 3 


4 (6.7) 


10 (76.9) 


4 (22.2) 




Atypical lymphocytosis 


8(13.3) 


7 (53.8) 


3 (16.7) 


Not observed 


Abnormal liver function 


7(11.7) 


11 (84.6) 


0(0) 




Abnormal renal function 


1 (1.7) 


2 (15.4) 


0(0) 





SJS/TEN, Stevens-Johnson syndrome/toxic epidermal necrolysis; HSS, hypersensitivity syndrome; MPE, maculopapular exanthema. 
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Fig. 1 

HLA-A HLA-C HLA-B HLA-DRB1 




Association screen of 20 short tandem repeat polymorphism (STRP) markers in the major histocompatibility complex (MHC) region with 
carbamazepine (CBZ)-induced cutaneous adverse drug reactions. On the jr-axis, 20 STRP markers in the MHC region are ordered by their physical 
positions (29.9-33.9 Mb) on chromosome 6p21.3. On they-axis, the -log 10 P values were calculated by comparison of the allele frequencies 
between the patients and tolerant group using the Cochran-Armitage exact trend test. Genotyping data of four groups of patients are presented: 
diamond symbols: CBZ-induced Stevens-Johnson syndrome (SJS)/toxic epidermal necrolysis (TEN); square symbols: CBZ-induced maculopapular 
eruption (MPE)/hypersensitivity syndrome (HSS); triangle symbols: CBZ-induced HSS; and circle symbols: CBZ-induced MPE. 



tolerant controls [OR = 254.6; 95% confidence interval 
(CI) = 70.4-901.9] for TT/TG genotype versus GG 
genotype) (Fig. 3). To verify the results of our association 
screen, we also genotyped 93 healthy subjects randomly 
selected from the general population [30], A strong 
association between SNPs and CBZ-induced SJS/TEN 
was replicated in the MHC region, especially rs3130690, 
nearby the HLA-B locus (data not shown). 

We further genotyped the individual HLA-A, B, C, DRB1 
and MICA alleles. The HLA-B*1502 allele was present in 
98.3% (59/60) of CBZ-SJS/TEN patients, whereas only 
4.2% (6/144) of the tolerant controls were positive for the 
allele (/» c =1.6x lO" 41 , OR = 1357; 95% CI = 193.4- 
8838.3) (Table 3 and Fig. 3). The only CBZ-SJS patient 
who did not have the B*1502 allele, had instead, another 
HLA-B 15 allele: HLA-B*1558. In addition to the positive 
association, we also found that HLA-B*4001 was negatively 
associated with CBZ-SJS/TEN (P c = 2.6 x 10 -4 , OR = 
0.16; 95% CI = 0.1-0.4) (Table 3). The HLA-Cw*0801 and 
MICA*019 alleles flanking the HLA-B* 1502 showed strong 
linkage disequilibrium with HLA-B* 1502, and were present 
in 93.3% (56/60) and 95% (57/60) of SJS/TEN patients, 
respectively. An extended B*1502 haplotype formed by 
polymorphic alleles (A*1101-Cw*0801-HLABC-CA*119- 
rs3130690T-B*1502-MICA*019-DRBl*1202) had a strong 
association with CBZ-induced SJS/TEN. The recombinant 
map of Cw*0801-HLABC-CA*119-rs3130690T-B*1502- 
MICA*019 defined the susceptible region within 86 kb 
(i.e. between the T allele of rs3130690 and MICA*019) 



flanking the B*1502 gene in the 4Mb MHC region 
(Fig. 3). Within this 86-kb region, HLA-B is the only 
known gene. Taken together, the data suggested that one 
or more alleles in the vicinity of the HLA-B locus, 
particularly B*1502 itself, participate in the pathogenesis of 
CBZ-induced SJS/TEN. 

Association screen for candidate gene SNPs with 
CBZ-MPE/HSS 

When MPE and HSS were grouped together as non- 
bullous cADRs, no STRP markers showed association, 
and only six SNPs had /'-values less than 0.01 (Figs 1 and 
2b). The three most significant SNPs are listed in 
Table 2, including rsl264511 located near HLA-E, 
rsl 042389 located on the 3'-untranslation region of 
CYP2B6, and rs2894342 located in the promoter region 
of the motilin (MLN) gene (Table 2). However, the 
P-values became non-significant after correcting for 
multiple testing (278 SNP assays) was performed. 

In HLA-typing, HLA-A*3101 showed an association with 
MPE/HSS (/>. = 0.0021) (Table 3). HLA-A*3101 was 
present in 25.8% (8/31) of patients with MPE/HSS, but 
only in 2.8% (4/144) of tolerant controls (OR =12.17; 
95% CI = 3.6-41.2). 

Association between MLN polymorphisms and HSS 

When genetic analysis was performed separately for HSS 
and MPE, one STRP marker, D6S2414 (physical map: 
32.9 Mb, chromosome 6) showed an association with 
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CBZ-induced HSS (/» = 6.25 x 10 -5 , Fig. 1). In addition, 
several SNPs located in the MHC class II terminal region, 
around the motilin gene (physical map: 33.8 Mb), showed 
associations with HSS (Table 2 and Fig. 2c). The most 
significant SNP rs2894342 was located in the promoter of 
the MLN gene, and its A allele had an increased risk for 
GBZ-HSS (P c = 0.0064, OR =7.11; 95% CI = 3.1-16.5) 
(Table 2). The next most significant SNP was a 
nonsynonymous SNP (rs2075800 of HSPA1L, heat shock 
70 kDa protein 1 -like) involving an A to G transition 
which leads to an amino acid change of lysine to 
glutamate at residue 602. The third most significant 
SNP was rs2395402, which was located in the intron of 
the LEM domain containing 2 (LEMD2) gene, 21 kb 
telomeric to rs2894342 of the MLN gene (Table 2). 
However, the /^-values for the latter two SNPs became 
non-significant after Bonferroni's correction. 

None of the HLA-A, B, Cand DRB1 alleles had significant 
associations with CBZ-induced HSS (Table 3). 



n HLA-A*3101 and MPE 

STRP or SNP markers in the HLA-B or MLN region on 
chromosome 6p21.3 showed no association with MPE 
(Figs 1 and 2d). Instead, two SNPs near the HLA-E 
region (physical position 30.5 Mb) and one SNP in the 
5'-untranslation region of the leukocyte specific transcript 1 
(LST1) gene showed associations with MPE 
(Fig. 2d and Table 2, P = 0.00078-0.0075). The two 
SNPs in the HLA-E region were: rsl 26451 1, an intergenic 
SNP located near HLA-E and rsl059510 located in HLA- 
E gene (Table 2). HLA-genotyping further revealed that 
the HLA-A*3101 allele, approximately 530 kb telomeric to 
HLA-E, was present in 33.3% (6/18) of MPE patients but 
only in 2.8% (4/144) in the tolerant control group 
(P c = Z.2x 10"\ OR =17.5; 95% CI = 4.6-66.5) 
(Table 3). 

Association screen for SNPs in drug metabolizing 
enzymes 

Genotyping of 88 informative SNPs in genes that are 
potentially involved in metabolic pathways for CBZ, 
including CYP3A4, CYP1A2, CYP2C8, CYP2B6, CYP2G9 
and EPHX1, in various types of CBZ-induced cADRs, 
showed that only one SNP, rsl 042389, in the 
3'-untranslation region of CYP2B6, had a weak association 
with CBZ-MPE/HSS (/> = 0.0016, Table 2 and Fig. 2b). 
The lvalue became non-significant after Bonferroni's 
correction, suggesting that the genetic polymorphisms of 
genes encoding drug metabolizing enzymes have no 
significant association with CBZ-induced cADRs, regard- 
less of the type of cADRs. 

Discussion 

This study confirmed our previous observation that HLA- 
B*1502 is strongly associated with CBZ-SJS/TEN in Han 
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Fig. 2 




Screening of 278 candidate single nucleotide polymorphisms (SNPs) for association with (a) carbamazepine (CBZ)-induced Stevens-Johnson 
syndrome/toxic epidermal necrolysis, (b) CBZ-induced maculopapular eruption (MPE)/hypersensitivity syndrome (HSS), (c) CBZ-induced HSS and 
(d) CBZ-induced MPE. On the A-axis, 278 SNPs are ordered by their chromosomal positions, including 58 SNPs of epoxide hydrolase 1 (EPHX1 ) on 
chromosome 1 , 1 90 SNPs in the major histocompatibility complex (MHC) region on chromosome 6, and 30 SNPs of CYP3A4, 2C8, 2C9, 1 A2 and 
2B6. On the y-axis, the -log 10 P-values were calculated by comparison of the allele (diamond symbols) or genotype (triangle symbols) frequencies 
between the patients and tolerant group using the Cochran-Armitage exact trend test. 



Chinese [18]. We increased the number of study patients 
from the original 44 to 60 and included patients from a 
wide geographical distribution (Taiwan, Hong Kong, 
China and the USA). However, all patients enrolled were 
Chinese or Chinese descendants. The fact that the HLA- 
B*1502 allele was present in low frequency in Caucasians 
may explain the low incidence of CBZ-SJS/TEN in 
Caucasians [31,32]. 

It is interesting to note that the only CBZ-SJS/TEN 
patient who did not have the HLA-B*1502 gene had 



HLA-B*1558, another HLA-B15 variant. This variant 
allele was present in low frequency in our population 
(approximately 0.9%) and was not detected in any of the 
144 tolerant patients. The detection of B*1558 in one of 
our patients may imply that B*1558 shares a similar 
structural feature with B*1502 for triggering the immune 
reaction of SJS caused by CBZ. 

Fine recombinant genetic mapping on the extended HLA- 
B*1502 haplotype (A*1101-Cw*0801-HLABC-CA*119-rs 
3130690TB*1502-MICA*019-DRB1*1202) further narrowed 



Copyright © Lippincott Williams & Wilkins. Unauthorized reproduction of this article is prohibited. 



304 Pharmacogenetics and Genomics 2006, Vol 1 6 No 4 



STRP marker 
HLABC-CA 



_J^JI325kb_^2^J_ 



rs31 30690 



Centromere 
HLA-DRB1 





A*1101 •■ 


Cw*oaoi 


HLABC 

-CA'119 


rs31 30690 B > 1502 

T allele 


MICA*019 -■ 


C 


SJS/TEN 


70% 


93.3% 


90% 


95% 


98.3% 


95% 


68.3% 


Tolerant 


46.5% 


13.9% 


4.2% 


6.9% 


4.2% 


1 3.2% 


16% 


P value 


0.003 


4x1 0~ 28 


2x10" 35 


4X10-36 


4x10" 43 


5x10- 30 


7x10" 




N.S 


8x10-" 


4x10" 34 


1x10- 35 


2x10" 41 


1x10" 28 


2x10" 


Odds ratio 


2.7 


86.8 


207.0 


254.6 


1357 


125 


11.4 


(95% C.I.) 


(1.4-5.1) 


(29.3-254.6) 


(65.0-659.1) 


(70.4-901.9 


(193.4-8838 


3) (37.3-412.4) 


(5.64- 



Susceptible region 

Susceptible region for carbamazepine (CBZ)-induced Stevens-Johnson syndrome (SJS)/toxic epidermal necrolysis (TEN) on the extended HLA- 
B*1502 haplotype. Schematic map of 4Mb major histocompatibility complex region on chromosome 6p21.3 shows the relative positions of HLA A, 
B, C, DRB1 and MICA genes, STRP marker (HLABC-CA), and SNP marker (rs31 30690). The percentage of presence of the specific allele in CBZ- 
SJS/TEN patients and tolerant-controls is shown. The susceptible region should be within 86 kb (i.e. between T allele of rs31 30690 and 
MICA*019). Within this region, HLA-B*1502 is the only known gene and shows the strongest association with CBZ-induced SJS/TEN. NS, Not 
significant (P>0.05). 



Table 3 Association of HLA alleles with carbamazepine (CBZ)-induced cutaneous adverse reactions 



SJS^EN (n=60) 



MPE/HSS (n=31) 



•4001 
HLA-Cw 
•0102 
•0801 



42 [NS; 2.68 (1.4-5.1) 
5 [0.026; 0.23 (0.1-0.6 
1 [NS; 0.59 (0.1-4.1) 

59 [1.6 xlO -4 '; 135^ 
(193,4-8838.3)] 
6 [2.6 x 10" 4 ; 0.16 (0.1-' 



12 [NS; 0.73 (0.3-1.6)] 

9 [NS; 1.03 (0.5-2.4)] 

8 [0.0021; 12.17 (3.6-41.2)1 

1 [NS; 0.77 (0.1-5.1)] 

10 [NS; 0.69 (0.3-1.5)] 



6 [NS; 0.99 (0.3-3)] 
3 [NS; 0.75 (0.2-2.7)] 

2 [NS; 6.36 (1.2-33.9)] 

0 [NS;0.79 (0.1-8.8)] 

3 [NS; 0.43 (0.1-1.5)] 



6 [NS; 0.5I 

6 [NS; 1.21 
6 [2.2 x 10 -3 ; 

1 [NS; 1.3! 

7 [NS; 0.9! 

9 INS; 1.6' 

2 [NS; 0.7! 



(0.2-9.3)] 
(0.4-2.4)] 



Data are genotype, n of positive subjects [P c value; odds ratio (95% confidence interval)] . The association of HLA-alleles was examined by Fisher's exa 
values were adjusted by using Bonferroni's correction for multiple comparisons (17 for HLA-A, 40 for HLA-B, 19 for HLA-C and 30 for HLA-DRB1). N 
(P>0.05). SJS/TEN, Stevens-Johnson syndrome/toxic epidermal necrolysis; MPE, maculopapular exanthema; HSS, hypersensitivity syndrome. 



down the susceptibility region for CBZ-SJS/TEN to 
within 86 kb flanking the HLA-B gene. In this region, 
HLA-B is the only gene present. Thus, it is most likely 
that B*1502 itself is the susceptibility gene for CBZ-SJS. 
Furthermore, it has been proposed that the molecular 
mechanism causing SJS/TEN involves CD8 + cytotoxic 



T-cell-mediated cell death [16]. The HLA-B gene can 
elicit immune responses by presenting endogenous 
antigens to the cytotoxic Tcells and it has been proposed 
that the drug and/or its metabolites may bind to the 
peptides which are then presented by the HLA-B allele 
and recognized by the specific T cells [33,34). Our 
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cytotoxic assays also supported the view that HLA-E*1502 
is functionally involved in the cytotoxic response 
mediated by the activated T cells in CBZ-induced SJS 
(unpublished data). The recent findings of HLA-B*5101 
in abacavir hypersensitivity [35] and HLA-B*5801 in 
allopurinol-induced severe cADRs [26] further support 
the important role of the HLA-B gene in these serious 
drug-induced adverse conditions. 

SJS and TEN are cADRs characterized by bullous lesions. 

We observed a single patient who developed a bullous 
fixed drug eruption while on CBZ. Interestingly, this 
patient was also positive for HLA-B*1502 (data not 
shown). However, the present study clearly showed that 
non-bullous cADRs caused by CBZ are not linked to the 
HLA-B*1502 marker, suggesting that the genetic associa- 
tion of CBZ-cADRs is phenotype-specific. It should be 
pointed out that our previous study of allopurinol induced 
cADRs showed that both SJS/TEN and HSS were 
associated with a same HLA-B allele [26]. However, in 
this study, CBZ-HSS is not linked to the HLA-B marker 
as SJS/TEN. At the present time, we do not know why 
this discrepancy occurs. Other genetic or environmental 
factors may contribute and lead to various spectrum of 
adverse reactions caused by different drugs. 

Although there is no direct evidence that MPE and HSS 
are the same disease spectrum with differences in 
severity, HSS and MPE may be grouped together as 
non-bullous cADRs because of the similarity of the 
cutaneous manifestation [2,5,19]. However, the present 
study suggests that HSS and MPE are two disease 
entities because they are linked to different genetic 
markers: CBZ-MPE was associated with HLA-A*3101, 
and CBZ-HSS with the motilin gene polymorphisms in 
the MHC class II terminal region (Figs 1 and 2 and 
Tables 2 and 3). We did not observe the association of 
CBZ-HSS with the TNF2 allele (rsl800629; TNF-308A 
allele; physical position: 31.6 Mb; Fig. 2c) or the ancestral 
haplotype 8.1 on the MHC, as reported previously in a 
study of CBZ-induced HSS in Caucasians [22]. This 
could be due to different study populations or different 
criteria in delineating the clinical phenotypes. In the 
present study, the SJS/TEN patients outnumber the 
patients with MPE as a reflection of the recruitment 
process because we recruited patients primarily from 
inpatients who had more severe diseases. However, in 
Han Chinese, patients with CBZ-SJS/TEN were indeed 
more frequently seen than CBZ-HSS (data not shown). 
Because of the small numbers of the patients enrolled in 
this study, further studies with a larger sample size will be 
needed to confirm this initial observation. 

Naisbitt et al. [15] cloned T cells of patients with CBZ- 
MPE/HSS and observed that T cells may recognize CBZ 
depending on the presence of HLA class II (DR/DQ)- 



matched antigen-presenting cells. Interestingly, in our 
SNPand STRP screen and HLA typing data, we observed 
associations between CBZ-induced HSS and polymorph- 
isms markers located terminal to MHC class II region. 
However, the association of HLA-A*3101 with MPE 
appears in disagreement with the observation that 
CD4+ T cells are the major cell type found in the skin 
lesions of MPE [15]. This could suggest that the gene 
associated with MPE lies in the vicinity of the HLA-A 
locus which is in linkage disequilibrium with HLA- 
A*3W1. 

In conclusion, our data show that genetic susceptibility to 
CBZ-induced cADRs is phenotype-specific. The suscept- 
ibility gene for CBZ-induced SJS/TEN lies within an 
86 kb region flanking the HLA-B locus and HLA-B* 1502 
itself could be directly involved in the pathogenesis of 
CBZ-induced SJS/TEN. The tight association of HLA- 
B*1502 and CBZ-SJS/TEN provides a plausible basis for 
the development of such a test to identify individuals at 
risk for this potentially life-threatening condition caused 
by CBZ in Han Chinese, as well as for a further increased 
understanding of the pathogenesis of the clinical 
syndrome. 
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